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Abstract. We propose a numerical method for finding solitary wave solutions of 
generalized Korteweg-de Vries equations by solving the nonlinear eigenvalue problem 
on an unbounded domain. The artificial boundary conditions are obtained to 
make the domain finite. We specially discuss the soliton solutions of the K(m,n) 
equation and KdV- K(m, n) equation. Furthermore for the mixed models of linear and 
nonlinear dispersion, the collision behaviors of soliton-soliton and soliton-antisoliton 
are observed. 
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1. Introduction 

In this paper, we propose a numerical method for finding solitary wave solutions, which 
we call numerical solitons, of generalized Korteweg-de Vries (GKdV) equations. The 
method is simple and effective to find numerical solitons and can thus be used to study 
the nonlinear mechanism of nonlinear evolution equations, especially in equations for 
which it is difficult to obtain soliton solutions by analytical tools. A sequence of GKdV 
equations which are not necessarily integrable are shown to admit soliton solutions by 
our method. In general, the GKdV equations are of the form 

u t + G(u, u x ) x + H (u) xxx = 0, (1) 

where G and H are given functions. The K(m,n) equation with m,n > 1, 

ut + {u m ) x + {u n ) xxx = 0, (2) 

is a special case of the GKdV equation (JTJ). 

When m = 2, n = 1, it is the first important model equation, the Korteweg-de 
Vries (KdV) equation, 

u t + {u 2 ) x + u xxx = 0, (3) 

which is integrable and admits sech 2 soliton solutions. For the KdV equation (JHJ), 
the linear dispersion and the nonlinear convection work against each other and exactly 
balance to result in a stable solution. Since Zabusky and Kruskal [1] numerically studied 
the elastic collisions of KdV solitons, the research of solitons has attracted many workers' 
attentions and interests. 

For n > 1, the K(m,n) equation(j2J) is fully nonlinear [2]. Their soliton solutions 
have compact support and called compactons. For m = n, they assume a very simple 
form [2], namely 

r 2\n 9r n — 1, , ... i x . nn 

u = {^i cos [-^(s-**)]}^ for \ x ~ Xt \^—[' ( 4 ) 

and zero otherwise, where A is a constant which represents the velocity of the travelling 
wave. The K(m,n) equation is not integrable [2,3]; however, the collisions of the 
compactons are almost elastic by numerical experiments [2]. This suggests that the 
elastic collisions are shared by many nonlinear equations, even though they are not 
integrable. Hence, it seems of great interest to observe more equations to further study 
their nonlinear mechanisms. 

However, the analytic ways are limited for many equations, due to the complexity of 
the nonlinear properties. The aim of the proposed numerical method is to obtain soliton 
solutions for these equations. Our approach is to transform the nonlinear eigenvalue 
problems on the unbounded domain deduced from the travelling wave solutions of the 
equations into bounded problems by imposing an artificial boundary condition. These 
problems can then be solved numerically in an efficient way. 

The rest of the paper is organized as follows. In section 2, we discuss the numerical 
method. In section 3, some application examples are proposed by using the method. 
Conclusions are given in section 4. 
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2. The method 

2.1. K(m,n) equations 

Although the method is suitable for wider classes of nonlinear evolution equations, here 
we introduce the idea by considering the solitary wave solutions of K(m, n) equation 
(J2J) for m, n > 1. Suppose that 

u(x,t) = f(0 = f(x-\ 1 t) (5) 

is a travelling wave solution of the K(m, n) equation. It is easy to see that Ai and /(£) 
satisfy the nonlinear eigenvalue problem 

- Ax/ + f m + (/»)« = o, (6) 

/(0 ^ as £ ^ ±oo. (7) 

In general, the spectrum is continuous, and if /(£) is an eigenfunction corresponding 
to Ai, then so is /(£ — £o) f° r an Y £o £ -R- Hence, it is very difficult to solve (6), (7) 
numerically 

If, however, we prescribe the amplitude of / and the location of its maximum, for 
example by setting 

/(0) = 1, (8) 

f(0) = 0, (9) 

then (6)-(9) represents a nonlinear boundary-value problem on (— oo,0] and [0,oo). 
Moreover, if Ai and /(£) is a solution of (6)-(9), then 

u(x 1 t) = Af(B(A)(x-\(A)t)), (10) 

is a solitary wave solution of equation K(m,n), with amplitude A > 0, and functions 
B(A) and X(A) to be determined. Let rj = B(A)(x - X(A)t). Substituting (JTUJ) into 
equation 0, and using equation (0), we obtain 

-AX(A)f + + A n B(A) 2 (f n ) m 

= (-AX(A) + A m X,)f + (A n B(A) 2 - A m )(f n ) m . 



(11) 



We set 



and 



m — n 

B(A)=A—, (12) 



X(A) = Aiv4 m . (13) 

We know that u(x,t) is a solitary wave solution of K(m,n) with amplitude A and 
velocity XiA m . Therefore, we only need to find Ai and /(£), an d this allows us to find 
the other solitary wave solutions of K(m,n). 

For the numerical solution of the problem (6)- (9) on the unbounded domain 
(— oo, +oo), we only need to consider the problem on the domain [0, oo), for the reasons 
of symmetry. We introduce an artificial boundary at the point £ = b > 0, where b is 
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large enough. We use the boundary condition f(b) = to reduce the problem to a 
problem on the bounded interval [0, b] for m > 1, n > 1: 

Find Ai G R and /(£), such that 

- Ai/ + / m + (/")« = 0, < e < 6 (14) 
/(0) = 1, (15) 
/'(0) = 0, (16) 
f(b) = 0. (17) 

We solve the problem numerically by a finite difference method. Let h = b/N with 
a positive integer N, ^ = ih, i = 0, 1, • • • , N. Then the interval [0, b] is divided into iV 
subintervals by 

= Co < 6 < • • • < &v = 6. 
Using the second-order central difference, we obtain 

« o in I in 

-A 1 / J + /r + Zm ^ — = 0, i = 0,l,..-,7V-l. (18) 

Due to the boundary conditions we have 

/-i = /i, /o = l, /at = 0. (19) 

The problem (18) with boundary conditions (19) is a nonlinear algebraic system which 
can be solved by the Newton iteration method. 

2.2. MKdV-K(m,n) equation 

The second model we consider here is the mKdV-K(m, n) equation 

u t + (u m ) x + (u + u n ) xxx = 0. (20) 

We also assume the travelling wave solution 

u(x,t) = f(£) = f(x-\t). (21) 

Then we have 

- A/< + (f m )i + (/ + Dm = 0, (22) 
which can be integrated once to yield 

- A/ + f m + (/ + = C, (23) 

where C is an arbitrary constant. Since we are interested in a localized solution, we 
know that: 

(i) . /(£) — > 0, as ^ — > ±oo; therefore we have C = in equation (23); 

(ii) . /(0) = A is an extremum, i.e., /'(0) = 0. 
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Similarly, we obtain the following nonlinear problem on [0, oo): 

- A/ + f m + (/ + D« = 0, (24) 

/(0 ^ as £ ^ +oo, (25) 

7(0) = A, (26) 

/'(0) = 0. (27) 

We introduce the artificial boundary point £ = 6 > 0, with 6 large enough, and m,n > 1. 
On the interval [6, oo), the equation (24) can be linearized to 

-A/ + /« = 0, £e[6,oo), (28) 



which has two linearly independent solutions, 

and e~^. 

By the condition (25), /(£) = Ce~^ with A > 0, and 

/'(0 = -Vxce-s* = -Va/(0- 

Then at £ = 6 we have the artificial boundary condition 



/'(&) = -VA/(6). (29) 
Thus we have reduced the problem to a problem on the bounded interval [0, b]: 

- A/ + /™ + (/ + /")€€ = 0> 0<£<6 (30) 

7(0) = A (31) 

/'(0) = 0, (32) 

f(b) = -V\f(b), (33) 
which can again be solved by the finite difference method. 



3. Numerical examples 

Example 1. Use the proposed method in Section II, numerical solitary wave solutions 
of the GKdV equation can be obtained. We first consider the K{2, 2) and K(3, 3) 
equations to test the efficiency of the method. 




Fig. 1. Numerical solitons. Solid line is the exact solution, (a) K(2, 2); (b) K(3, 3). 
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In this case, the Newton iteration maybe does not converge or the solution shows 
blowup phenomena with the grid refinement. To avoid these phenomena, a low pass 
filter has been used in iterate steps. Fig. 1 (a) and (b) show the numerical results with 
h = 0.25. We see the numerical results are well agreed with the exact solutions. 

Example 2. We show the numerical soliton to the KdV-_fT(2, 2) equation 

ut + (2m 2 ), + (u + u 2 ) xxx = 0, (34) 

which is a combination of the KdV and K(2, 2) equations. Unlike the KdV (j3J) and 
K(m, n) (J2J equations, it is a mixed dispersive model. Both linear and nonlinear 
dispersions are often presented in many physical applications, for example the Carmassa- 
Holm equation [4] in water waves. It is interesting to solve them to study the interaction 
mechanism of linear and nonlinear structures [5] . 

Fig. 2(a) illustrates the shapes of three numerical solitons with amplitudes A = 2, 
1.5 and 1, respectively. We also perform the numerical experiments of collisions with 
the three solitons which are located at the three isolated centers x = —60, -30 and 0. 
Here, the fully implicit Crank-Nicolson scheme is used in the domain [—80, 80] with 
periodic boundary conditions. Fig. 2(b-d) shows the evolution of the three numerical 
solitons. Extensive tests all indicate that the collisions of the KdV-K(2, 2) equation 
(jnij) are almost elastic. 



(a) 



(b) 





(d) 




Fig. 2. The evolution of three numerical solitons. 



Similarly to the behavior of the K(2, 2) compactons, the collisions of KdV-f^(2, 2) 
numerical solitons also create a small ripple (see Fig. 3). The amplitudes of the ripple 
are less than 3% of the amplitudes of all the solitons. However, their amplitudes do not 
almost decrease with the mesh is refined. The result illustrates that the creation of the 
collision ripple probably does not arise from the numerical accuracy but from its internal 
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mechanism of the equation. Though this possibility appears reasonable, it needs to be 
confirmed by more proofs since now it has only been deduced from numerical results. 




-10 10 20 30 40 50 60 -10 10 20 30 40 50 60 



(a) * (b) 

Fig. 3. The ripples created by the collision of solitons in different grids. 

Example 3. The model equation we considered is called the mKdV-K(3, 3) 
equation, 

u t + (2u 3 ) x + (u + u 3 ) xxx = 0, (35) 

which is a combination of the modified KdV and K(3, 3) equations. Since all terms 
of the mKdV-K(3, 3) equation (|3*3|) are odd, its soliton solutions are in pairs by the 
relationship (u, —u). They propagate in the same direction with the same speed. The 
solitons with negative signs are normally called antisolitons. 

The K (2, 2) equation also admits anticompactons which propagate in the opposite 
direction. However, it is difficult to numerically simulate compacton-anticompacton 
collisions because of instability. It is not clear whether or not the reason of this instability 
is due to numerical errors [2]. Thus, in order to better understand the reason, it is 
meaningful to first study the behavior of soliton-antisoliton collisions. 

Fig. 4 shows the evolution of numerical soliton and antisoliton with amplitudes 2 
and 1, respectively. Here we take the spatial mesh size h = 0.025 in [-40,40] and periodic 
boundary conditions. We can see that their collisions are elastic if the tail oscillations 
originate from numerical in accuracies. In fact, extensive numerical experiments indicate 
that the oscillations decrease when the mesh is refined. This illustrates that the creation 
of oscillations are mainly caused by the numerical errors. 




(a) * (b) 
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(c) * (d) 

Fig. 4. The evolution of soliton-antisoliton collisions. 

These oscillations are high-frequency waves which are caused by every collision 
because of the numerical errors. We see that the collisions of oscillations and solitons 
are also elastic (see Fig. 4(c) (d)). However, it is difficult to handle them unless we use 
a very fine grid size. This is also the reason why compacton-anticompacton collisions of 
the K(m, n) equation are difficult to be simulated. Based on the above experiments, it 
is conjectured that the K(m, n) compacton-anticompacton collisions are also elastic. 

We know the Benjamin-Bona-Mahoney (BBM) equation 

iH + {u 2 ) x - u xxt = (36) 

as a good example of a nonintegrable model that has also soliton and antisoliton 
solutions. Its soliton-soliton collisions are elastic, but numerical tests illustrate that 
its soliton-antisoliton collisions are inelastic [6]. Therefore, it is interesting to study 
what governs the behavior of solitons collisions, but we are still a long way from a good 
understanding of the wave mechanics [7]. 

4. Conclusion 

In summary, we have reported a method to find numerical solitons of the GKdV 
equations which are not necessarily integrable. Using this method, we have given 
a basic study of soliton-soliton collisions of the KdV-fT(2, 2) equation and soliton- 
antisoliton collisions of the mKdV-i^(3, 3) equation. This method and the numerical 
experiments will contribute to a future, more extensive study of the theoretical aspects 
of the nonlinear dispersive mechanism of evolution equations. 
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